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Abstract 

We represent a theory of polymer gelation as an analogue of liquid-glass transition 
in which elastic fields of stress and strain shear components appear spontaneously as 
a consequence of the cross-linking of macromolecules. This circumstance is explained 
on the basis of obvious combinatoric arguments as well as a synergetic Lorenz sys- 
tem, where the strain acts as an order parameter, a conjugate field is reduced to 
the elastic stress, and the number of cross-links is a control parameter. Both the 
combinatoric and synergetic approaches show that an anomalous slow dependence 
of the shear modulus on the number of cross-links is obtained. 
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1 INTRODUCTION 

Within the phenomenological framework the basic distinction between a liquid 
and a glass consists in character of the relaxation law of shear components of 
the elastic stress: if in an ideal glass they are kept infinitely long, in a liquid 
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such relaxation proceeds in the finite time t = rj/G^ where rj is the dynamical 
shear viscosity and G is the shear modulus [1]. In naive manner it is possible 
to assume, that a glass transition is caused purely by kinetic effect of a liquid 
freezing, for which viscosity r] gets infinite value for a finite value of the shear 
modulus G [2]. However in the course of usual second-order phase transition, 
where infinite increase of the r at critical point is also observed the situation 
is reverse. Really, proceeding from the viscoelastic liquid to a general case, 
one has t = x/li where x is the generalized susceptibility and 7 is the kinetic 
coefficient (in the case under consideration they are reduced to the quantities 
G^^, r]^^, respectively) [3]. The infinite increase of a susceptibility x occurs 
and a kinetic coefficient 7 does not manifest any peculiarity at the phase 
transition. In our case, this is equivalent to the situation when the shear 
modulus G is approaching zero value at a fixed viscosity rj. This situation 
corresponds to a viscoelastic transition [4] . 

Usually under a glass transition, the transition into a glassy state which occurs 
during fast cooling of a liquid is understood. However, an analogous type of 
transition occurs during polymer network formation {polymer gelation) from 
monomers or linear macromolecules {sot) by means of a chemical reaction at 
constant temperature. This transition is referred to as the sol-gel transition. 
In such thermodynamic peculiarities are observed as, e.g., appearance 

of the shear modulus when the system reaches a critical point .[^ The first 
theoretical representations of such type transitions have been elaborated in 
classical works by Flory [5] - [7] . Modern description of the glassy state of the 
macro molecule networks (see [8] - [10]) is based on pioneering contributions 
to the theory of soft condensed matter: the Deam-Edwards theory [11] of a 
cross-linked tangled macromolecule, and the Edwards- Anderson theory [12] of 
spin glass. 

The aim of this paper is to elaborate the simplest theoretical scheme so that 
the transition occurring during polymer gelation can be explained on the base 
of obvious combinatoric arguments (Section 2) as well as within framework 
of a synergetic theory (Section 3). In the former case, the critical point and 
behaviour of the system behind the critical point are obtained using informa- 
tion about the cross-linking properties of the constituent parts of the system 
{monomers) and extent of the cross-linking process, only. In Section 3, the 
starting point is that a sol-gel transition is ensured by self-organization of the 
elastic fields of stress and strain shear components and the number of cross- 
links. Both the combinatoric approach and the synergetic one allow to obtain 
anomalous dependence of the shear modulus on the number of cross-links. 



Of course, when the temperature of the polymer network formed is decreased 
the kinetic glass transition in above sense occurs, as well. 
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2 Combinatoric approach 



The gelation process has been described first by Flory[5] and Stockmayer[13]. 
They used a laborious approach based on combinatoric considerations of the 
most probable composition of the system. A much more effective variant of 
this approach, exploiting theory of branching processes, has been developed 
later by Dobson and Gordon [14,15]. This approach allows also an analysis of 
structural details of the system behind the critical point, see, e.g. [16]. 

As the simplest case let us consider the gelation process in the system which 
at the beginning consists of a large number of monomers A/", wearing / func- 
tional terminal groups of independent reactivity. Monomers react mutually 
and irreversibly via their terminal groups and if / > 2, branched molecules of 
increasing size and complexity are formed progressively in the system. Extent 
of the reaction is described by the conversion of terminal groups a, which 
is defined as the ratio of the number of groups consumed by the reaction at 
given time and the starting number of groups. Eventually, when a critical con- 
version etc is achieved a molecule of macroscopic dimensions {gel) appears in 
the system [5] , [6] . Gel contains cycles of long sequences of linked monomers 
and, consequently, attains elastic properties. Flory [7] has shown that shear 
modulus of gel G is proportional to its "cycle rank" ^ defined as the number of 
"superfluous" links formed in gel which can be cut without breaking integrity 
of gel. In other words, cycle rank is the number of cuts needed for elimination 
of all cycles from gel. 

The key role in the approach is played by the extinction probability v, which 
is a probability that a link formed in the process has just a flnite continuation. 
It can be shown easily [14] that for the system considered here the extinction 
probability can be obtained as a root of the equation 



satisfying the condition < ^; < 1. Eq.(l) expresses fact that a given link has 
a finite continuation in a direction only if / — 1 functional groups remaining on 
the monomer connected by the link are either unreacted (with probability 1 — 
a) or reacted giving rise to links with finite continuation only (with probability 
av). Below the critical conversion, the only solution of Eq.(l) is v — 1, i.e., 
only molecules of finite size are formed in the system. However, behind the 
critical conversion monomers can be found either in sol or in gel: v < 1. As 
a measure of the "distance" from the critical conversion let us introduce the 
parameter e defined as 



V — {1 — a + av) 




e = 



a — 



(2) 
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and expand the extinction probability in a series: 

V = 1 + Ae + Ae^ + ... (3) 



Substituting Eqs.(2) and (3) into Eq.(l), expressions for the critical conversion 
and the parameters Aj, are obtained as 



1 



/-I' 

if -I? 



f-2 ' 
4 (2/-3)(/-lf 
3 (/-2)2 • 



(4) 
(5) 
(6) 



The expression for the cycle rank can be found in the following way. Obviously, 
to join J\f monomers into a cycle-free structure, J\f —1 links are necessary. By 
definition, the cycle rank ^ is the number of "superfluous" links formed in gel, 
i.e., the difference between the number Nq of all links formed in gel and the 
number Ag — 1 of links sufficient to join together A/q monomers of gel: 

^^Ng-(A/"g-1)^Ng-AAg (7) 



as Ng, a/" ^ 1. Correspondingly, the number of links formed in gel is the 
difference between the numbers of links formed in the total system and in sol, 
i.e., 

Ng = ^A^«/(1-^^) (8) 



as any link in sol has to have finite continuations in two directions. On the 
other hand, the number Aq of monomers incorporated in gel is the difference 
between the total number of monomers and the number of monomers in sol 
which is made of monomers with links of only finite continuation: 

j\fQ^j\f-j\f(^l-a + avy. (9) 



By virtue of Eqs.(7) - (9), one gets for cycle rank of the system considered 
i. ^ (1 _ « + avY + ^af{l - v') - 1. (10) 
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Finally, substituting Eqs.(3),(5) and (6) into the formula (10), one gets the 
necessary expansion: 

^ 2 (/--l)^ 3 4 

U 3 (/-2)2-^ ^ ^ >- ^ ' 

Respectively, the weight fraction of gel wq = Mq/H is determined by Eq.(9) 
to read 

wo^l^-l^e + Oie'). (12) 



So, if the critical exponent of the gel weight fraction is equal 1 as usual, this 
for the cycle rank is anomalous large being equal 3. It is worthwhile to note 
that combinatoric approach is based exclusively on the information about the 
functionality of the monomers and extent of the chemical reaction between 
the monomers in the system considered. 



3 Synergetic approach 



Now, let us consider the polymer network as a viscoelastic continuum matter 
that is characterized by the shear modulus G and the shear viscosity rj. Process 
of polymer gelation is determined by the number of the cross-links A^, which 
value is different from a stationary magnitude A^o at a time t. Therein, an 
elastic state of the polymer is defined by the shear component of the proper 
(internal) values of deformation e{t) and stress o'(t). The keypoint is that these 
values are not reduced to the external elastic deformation e -C 1 and stress 
(Te <^ G, in particular they can get large magnitudes £ ~ 1, cr ~ G. 

Our consideration of evolution of the elastic continuum with the internal struc- 
ture is stated on the phenomenological equations by Maxwell-Kelvin [1] 

^ = -- + -, (13) 
at T T] 

— = + g,eN. (14) 

at T„ 



Here wc introduce a macroscopic relaxation time r for the strain and a mi- 
croscopic one To- for the stress, as well as a constant g„ > Q oi the positive 
feedback between the deformation e and the number of cross-links N. Within 
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the microscopic interval t ^ t^, steady-state condition da /dt — in Eq.(14) 
leads to the Hooke law with the microscopic shear modulus 

= T^g^N (15) 



being determined by the number of cross-links N. Respectively, within a 
macroscopic interval t ':$> t Eq.(13) gives the magnitude G = tj/t that is 
characteristic for the usual modulus of the viscoelastic matter. Lastly, a varia- 
tion rate dN/ dt of the internal degree of freedom is supposed to be determined 
by the equation 

dN Nq-N 

— = gN(T£ (16) 

at Tn 



where rjv is a mesoscopic relaxation time, 5^Jv > is constant of negative 
feedback between the deformation e and the stress a. Within a mesoscopic 
interval tat "C t <C r, Eq.(16) determines a steady-state value 

N ^ No-TNgNcre (17) 



that is smaller than the magnitude A^o fixed by external conditions due to the 
fact that the elastic energy is proportional to the product ae. 

System of Eqs.(13), (14) and (16) is known in synergetics [17] as the Lorenz 
system where the deformation e, the stress a and the number of cross-links N 
play roles of an order parameter, a conjugate field and a control parameter, 
respectively. It is very important for following considerations that the relation 
between micro-, meso- and macroscopic values of the relaxation times 

Ta,r7V<T (18) 



is satisfied. Due to this condition the evolution of the quantities a, N turns 
out to be subordinated to the long-time variation of £. A peculiarity of the 
Lorenz system consists in linear character of the equation (13) for the order 
parameter e and in non-linearity of equations (14), (16) for the conjugate field 
a and the control parameter N. The negative nature of non-linearity in Eq.(16) 
means a decrease of the number of cross-links. Evidently, this fact reflects 
Le Chatelier principle. A non- linear term in Eq.(14) for a field a describes the 
positive feedback causing the system self-organization. 

Expressions (13), (14) and (16) form the complete system of equations de- 
termining the polymer cross-linking behaviour. Because of a slow evolution, 
the order parameter e{t) subordinates variations of quantities cr{t), N{t), so 
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that one can take da/dt = dN/dt = within the framework of the adiabatic 
approximation [17]. Then A^, a are expressed in terms of e by the equations: 

^ ^ 1 I ^2/ 2 ' ^rn^ = ^<^^N9a9N; (19) 

^ = , ^0 / 2 ' = r^^^aATo. (20) 



In accordance with Eq.(19) the number of cross-hnks decreases monoto- 
neously with increase of the strain e from the value iVo at e = to Ai"o/2 
at e = In Eq.(20) the elastic stress in terms of the strain has the lin- 

ear form of the Hooke law at e <^ with the effective shear modulus Go- 
Then, at e = Em the function a{e) has a maximum and at e > it decreases 
which is physically meaningless. Thus, the constant defined by the second 
equation (19), has the meaning of the maximum achievable strain. 

Substituting Eq.(20) in Eq.(13) we find the equation describing evolution of 
a system in the course of the sol-gel transition: 



de 
'di 



dE 



7 



tT No 



(21) 



where constant 7 plays a role of the kinetic coefficient. Behaviour of the system 
under consideration is determined by the dependence E{e) of the elastic energy 
on the strain: 



E 



T Nn 




(22) 



where the characteristic value of the number of cross-links is introduced 

V 



(23) 



At A^o < dependence (22) is monotoneously increasing with a minimum at 
the point e = 0. It means that in the stationary state (^=0) the elastic strain 

^ Obviously this decrease is caused by the negative feedback in Eq.(16), that is 
reflection of Le Chatelier principle for analyzes problem. Actually, a liquid self- 
organization, resulting in a sol-gel transition, is caused by the positive feedback 
between the strain e and the number of cross- links in Eq.(14). Hence, the increase 
of the number of cross-links should intensify the self-organization effect. However, 
according to Eq.(19) system behaves in such way that the consequence of self- 
organization, i.e., growth of the elastic strain, leads to decrease of its origin - the 
number of cross-links. 
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is not spontaneous. Thus, a liquid state is realized, in which the strain caused 
by the external stress relaxes during the time 

Tef = r{l-No/N,)-\ (24) 



The relaxation time increases infinitely when the number of cross-links Nq 
reaches the critical value Nc and at Nq > Nc the system undergoes a sol-gel 
transition. In the gel state the multiplier 1/2 appears in the dependence (24), 
and the minimum of the elastic energy corresponds to the elastic strain 



r2 



A'o N, 1 + e 



where we introduce the distance from the critical value Nr 



being equivalent to the definition (2). Inserting Eq.(24) into the dependence 
(22), we obtain the elastic energy of the steady-state: 

T/V 

Eo^Eieo) = - ^e' + Oie'). (27) 



As would be expected, this energy is proportional to the second power of the 
parameter (26) and is negative in nature (the latter means the energy benefit 
of the gel state in comparison with the liquid state). 

Taking into account that the glassy state is determined by density of the 
localized monomers, let us find now the shear modulus of the appeared gel 
state. It is principally important in our considerations that the gel state is 
determined by the value of the Deam-Edwards parameter of localization u 
[11] which is supposed to be proportional to the square of the proper strain 
Eo of the matter. Under the condition of the appearance of the elastic strains 
e, a generalized Deam-Edwards parameter uj{e) has to be considered which 
is related to the condition uj{e = 0) = a;. Then, expanding the function uj{e) 
into a series and keeping the first two terms only, one obtains 

u;(e) ~ u;(l + e^), on oc el, e < 1. (28) 



By virtue of the parity condition uj{e) = uj{—e), this expansion does not 
contain a linear term.Q Because of that the total strain e oc ^Juj{e), internal 

Otherwise, the usual squared dependence (31) will not be obtained. 
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one Eq oc y/cu and elastic one e are connected by the following relation: 

e^ = el{l + e^), e < 1. (29) 



The key point is that Eq.(29) supposes the additivity rule holds not for quan- 
tities e, Eq, e themselves, but for their squares. The physical reason for such 
a situation is that the system under consideration is random in character and 
described by a symmetrical distribution function. Therefore, all odd-power 
moments vanish identically and making use of the additivity rule (29) for 
variances follows. 

It is worthwhile to note the seeming contradiction between above relations 
for the localization parameter cu and the Deam- Edwards results [11]. Obvious 
reason consists in that the formers are obtained within framework of the mean- 
field theory, whereas the latters suppose fluctuation effects. According to [11] 
the cross-linking process is not sensible to strain e and corresponding local- 
ization parameter uu is proportional to the cross-link number N but not to the 
difference N — Nc, as in Eqs.(25), (26). In our opinion, this is caused by non- 
self-consistency of the approach [11] in sense that the stress field a is switched 
off. On the other hand, making use of the statistical scheme [11] arrives at the 
strain-dependence for the localization parameter uj due to appearance of the 
polymer network entanglement, whereas within the above phenomenological 
approach this dependence has to be postulated. 

A contribution to the elastic energy caused by the external strain is determined 
by equality 

AE{e) = \E{e{e))-E{eo)\. (30) 



With use of Eqs. (22), (27), (29) and expansion hi{l+x) ^ x-x^/2+x^/3, x <C 
1 the expression 
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AE{e) ^ No (31) 



is easily obtained where only the first non-zero term is kept. Comparing this 
relation with the usual expression for the elastic energy [1] 

AE{e) ^ V (32) 



where V is volume, we arrive at the final expression for the shear modulus of 
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the gel state of the polymer network: 



G 



< 1 



V_ 

iVo 



(33) 



e = 



The notable peculiarity of this result consists in that, in accordance with 
previous result (11), the shear modulus is proportional to the third power of 
the distance e. It is worthwhile to note that an expression of this kind can be 
obtained only within framework of the synergetic approach used, but not on 
the basis of the phase transition theory by Landau. 

Finally, comparing Eqs. (11) and (33), we obtain the relation between the 
micro- and macroscopic parameters of the gel under consideration in the prox- 
imity of the critical point (e <C 1): 
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